Stationary Josephson effect in a weak-link between nonunitary triplet superconductors 
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A stationary Josephson effect in a weak-link between misorientated nonunitary triplet super- 
conductors is investigated theoretically. The non-self-consistent quasiclassical Eilenberger equation 
for this system has been solved analytically. As an application of this analytical calculation, the 
current-phase diagrams are plotted for the junction between two nonunitary bipolar /—wave su- 
perconducting banks. A spontaneous current parallel to the interface between superconductors has 
been observed. Also, the effect of misorientation between crystals on the Josephson and spontaneous 
currents is studied. Such experimental investigations of the current-phase diagrams can be used to 
test the pairing symmetry in the above-mentioned superconductors. 
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I. INTRODUCTION 

In recent years, the triplet superconductivity has be- 
come one of the modern subjects for researchers in 
the field of superconductivity 0, 0- Particularly, 
the nonunitary spin triplet state in which Cooper pairs 
may carry a finite averaged intrinsic spin momentum 
has attracted much attention in the last decade 0, 
A triplet state in the momentum space k can be de- 
scribed by the order parameter A(k) = i(d(k) • d)o y 
in a 2x2 matrix form in which <jj are 2x2 Pauli matri- 
ces (<r = (a x , <7 y , cr z )). The three dimensional complex 
vector d(k) (gap vector) describes the triplet pairing 
state. In the nonunitary state, the product A(k)A(k)'l' = 
d(k) • d*(k) + i(d(k) X d*(k)) • a is not a multiple of the 
unit matrix. Thus in a non-unitary state the time rever- 
sal symmetry is necessarily broken spontaneously and a 
spontaneous moment m(k) = id(k) x d*(k) appears at 
each point k of the momentum space. In this case the 
macroscopically averaged moment < m(k) > integrated 
on the Fermi surface does not vanish. The value m(k) is 
related to the net spin average by Tr[A(k)^(7 : jA(k)]. It 
is clear that the total spin average over the Fermi surface 
can be nonzero. As an application, the nonunitary bipo- 
lar state of /—wave pairing symmetry has been consid- 
ered for the B— phase of superconductivity in the U Pt$ 
compound which has been created at low temperatures 
T and small values of the magnetic field -ff [H, |(| . 

In the present paper, the ballistic Josephson weak link 
via an interface between two superconducting bulks with 
different orientations of the crystallographic axes is in- 
vestigated. This type of weak link structure can be used 
for the demonstration of the pairing symmetry in the su- 
perconducting phase 7] . Consequently, we generalize the 
formalism of paper y| for the weak link between triplet 
superconducting bulks with a nonunitary order parame- 
ter. In the paper |§j], the Josephson effect in the point 
contact between unitary /—wave triplet superconductors 
has been studied. Also, the effect of misorientation on 



the charge transport has been investigated and a spon- 
taneous current tangential to the interface between the 
/—wave superconductors, has been observed. 

In this paper the nonunitary bipolar /—wave model 
of the order parameter is considered. It is shown that 
the current-phase diagrams are totally different from the 
current-phase diagrams of the junction between the uni- 
tary triplet ( axial and planar) /—wave superconductors 
[8j . Roughly speaking, these different characters can be 
used to distinguish between nonunitary bipolar /—wave 
superconductivity and the other types of superconductiv- 
ity. In the weak-link structure between the nonunitary 
/—wave superconductors, the spontaneous current paral- 
lel to the interface has been observed as a fingerprint for 
unconventional superconductivity and spontaneous time 
reversal symmetry breaking. The effect of misorienta- 
tion on the spontaneous and Josephson currents is in- 
vestigated. It is possible to find the value of the phase 
difference in which the Josephson current is zero but the 
spontaneous current, which is produced by the interface 
and is tangential to the interface, is present. In some con- 
figurations and at the zero phase difference, the Joseph- 
son current is not generally zero but has a finite value. 
This finite value corresponds to a spontaneous phase dif- 
ference which is related to the misorientation between 
the gap vectorsd. 

The arrangement of the rest of this paper is as follows. 
In Sec[n]we describe the configuration that we have in- 
vestigated. For a non-self-consistent model of the order 
parameter, the quasiclassical Eilenberger equations \H\ 
are solved and suitable Green functions have been ob- 
tained analytically. In Sec lIIII the formulas obtained for 
the Green functions have been used for the calculation of 
the current densities at the interface. An analysis of nu- 
merical results will be presented in Sec IIVI together with 
some conclusions in SecfVl 
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FIG. 1: Scheme of a flat interface between two superconduct- 
ing bulks which are misorientated as much as a. 



II. FORMALISM AND BASIC EQUATIONS 

We consider a model of a fiat interface y = between 
two misorientated nonunitary /—wave superconducting 
half-spaces (Fig.l) as a ballistic Josephson junction. In 
the quasiclassical ballistic approach, in order to calculate 
the current, we use "transport-like" quasiclassical Eilen- 
berger equations for the energy integrated Green 
functions g (vj?, r, e m ) 



Vf • Vg 



e m a 3 + iA,g =0 



(1) 



and the normalization condition gg = 1, where e m — 
irT(2m + 1) are discrete Matsubara energies m = 
0,1,2,..., T is the temperature, Vi? is the Fermi veloc- 
ity and <X3 = 03 (£> / in which &j (j = 1,2,3) are Pauli 
matrices. The Matsubara propagator g can be written in 
the form: 
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where the matrix structure of the off-diagonal self energy 
A in the Nambu space is 



d • <Tl<72 

icr 2 d* ■ a 



(3) 



The nonunitary states, for which d x d* ^ 0, are inves- 
tigated. Fundamentally, the gap vector (order parame- 
ter) d has to be determined numerically from the self- 
consistency equation U, while in the present paper, we 
use a non-self-consistent model for the gap vector which 
is much more suitable for analytical calculations So- 
lutions to Eq. Q must satisfy the conditions for the 
Green functions and the gap vector d in the bulks of the 
superconductors far from the interface as follow: 



1 

g= n~ 



e m (l 
i&2 [id* n 4 



- A„ • a) [id n - d n x A„] • 0102 
d* x A„] • cr -£<7 2 (1 + A„ • <t)<7 2 

(4) 



where 



id n x d; 



(5) 




2[(el + d n ■ d*)2 + (d„ x d*)2] 



el + d„ • d* + V(4 + dn-d;) 2 + (d n xd;) 2 



d(±x) = d_> .,. (T. v, )rxi» ( Ty 



(6) 
(7) 



where (f> is the external phase difference between the order 
parameters of the bulks and n = 1,2 label the left and 
right half spaces respectively. It is clear that poles of the 
Green function in the energy space are in 



Oi = 0. 



Consequently, 



(-£ 2 + d n .d;f + (d„xd;) 2 = o 



and 



E = ±Jd n ■ d* ± id n xd* n 



(8) 



(9) 



(10) 



in which E is the energy value of the poles. The Eq. 
l|T| has to be supplemented by the continuity condi- 
tions at the interface between superconductors. For all 
quasiparticle trajectories, the Green functions satisfy the 
boundary conditions both in the right and left bulks as 
well as at the interface. The system of equations Q 
and the self-consistency equation for the gap vector d 
T] can be solved only numerically. For unconventional 
superconductors such solution requires the information 
about the interaction between the electrons in the Cooper 
pairs and the nature of unconventional superconductiv- 
ity in novel compounds which in most cases are unknown. 
Also, it has been shown that the absolute value of a self- 
consistent order parameter is suppressed near the inter- 
face and at the distances of the order of the coherence 
length, while its dependence on the direction in the mo- 
mentum space almost remains unaltered [T^ . This sup- 
pression of the order parameter changes the amplitude 
value of the current, but does not influence the current- 
phase dependence drastically. For example, it has been 
verified in Refs. 10] for the junction between unconven- 
tional ci-wave superconductors, in Ref.^3 f° r the case 
of unitary "/-wave" superconductors and in Ref.|l3| for 
pinholes in 3 _ffe, that there is good qualitative agreement 
between self-consistent and non-self-consistent results for 
not very large angles of misorientation. It has also been 
observed that the results of the non-self-consistent model 
in |l4j are similar to experiment |l5j . Consequently, de- 
spite the fact that this solution cannot be applied di- 
rectly to a quantitative analysis of a real experiment, only 
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FIG. 2: Component of the current normal to the interface 
(Josephson current) versus the phase difference (f> for the junc- 
tion between nonunitary bipolar /—wave bulks , T/T c — 0.15, 
geometry (i) and different misorientations. Currents are given 
in units of j = §eN(0)v F A (0). 



a qualitative comparison of calculated and experimental 
current-phase relations is possible. In our calculations, a 
simple model of the constant order parameter up to the 
interface is considered and the pair breaking and scatter- 
ing on the interface are ignored. We believe that under 
these strong assumptions our results describe the real sit- 
uation qualitatively. In the framework of such a model, 
the analytical expressions for the current can be obtained 
for a certain form of the order parameter. 



III. ANALYTICAL RESULTS 

The solution of Eq. allows us to calculate the cur- 
rent densities. The expression for the current is: 



FIG. 3: Component of the current normal to the interface 
(Josephson current) versus the phase difference cj> for the junc- 
tion between nonunitary bipolar /—wave bulks , T/T c = 0.15, 
geometry (ii) and different misorientations. 

tic case the system of equations for functions g t and gi 
can be decomposed into independent blocks of equations. 
The set of equations which enables us to find the Green 
function g\ is: 

w F kV<?i = j(d-g 3 -d*-g 2 ); (12) 

U F kVg_ = -2(d xg 3 + d*xg 2 ); (13) 

i> F kVg 2 = -2e m g 2 + 2igid + d x g_; (14) 

■u F kVg 3 = 2e m g 3 - 2i gi d* + d* x g_; (15) 

where g_ = gi — g4. The Eqs. (|12|) - (|14fl can be solved 
by integrating over the ballistic trajectories of electrons 
in the right and left half-spaces. The general solution 
satisfying the boundary conditions at infinity is 



j (r) = 2ineTN (0) ^ (v Fjl (v F , r, e m )) 



(11) 



9i 



(n) 



a n e 



(16) 



where (...) stands for averaging over the directions of an 
electron momentum on the Fermi surface v F , and N (0) 
is the electron density of states at the Fermi level of en- 
ergy. We assume that the order parameter is constant 
in space and in each half-space it equals its value far 
from the interface in the left or right bulks. For such 
a model, the current-phase dependence of a Josephson 
junction can be calculated analytically. It enables us to 
analyze the main features of current-phase dependence 
for any model of the nonunitary order parameter. The 
Eilcnberger equations Q f° r Green functions g, which 
are supplemented by the condition of continuity of solu- 
tions across the interface, y = 0, and the boundary con- 
ditions at the bulks, are solved for a non-self-consistent 
model of the order parameter analytically. In the ballis- 
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(17) 



(n) *d n d„ x A n 2m„d„ + d n x C n _2sQ n t. 

g 2 — — ~ — ~ e 



(18) 



( „) _ »d t * + d T * x A n 2ia n d* n - d t * t x C n _ 2sQnt 
g 3 — ^ t ^ e 



(19) 

where t is the time of flight along the trajectory, sgn (t) = 
sgn (y) — s and n = sgn (v y ) . By matching the solutions 
(|16ll9j) at the interface (y = 0,t = 0), we find constants 
a n and C„. Indices n — 1,2 label the left and right half- 
spaces respectively. The function gi (0) = g^' (— 0) = 
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FIG. 4: The a;— component of the current tangential to the 
interface versus the phase difference <j> for the junction be- 
tween nonunitary bipolar /—wave superconducting bulks, 
T/T c — 0.15, geometry (i) and the different misorientations. 

g\ (+0) which is a diagonal term of the Green matrix 
and determines the current density at the interface, y = 
0, is as follows: 

v (d 2 ■ d 2 (yn 1 + ef - di ■ di(7?n 2 - ef + B) 
■ 9UUj [dafafix+ej+di OA -e)] a 

(20) 

where B = id x x d 2 ■ (A 1 +A 2 )(ryf2 2 _ e)(?7fli + e). We 
consider a rotation R only in the right superconductor 
(see Fig.l), i.e., d 2 (k) = .Rdi(.R _1 k); k is the unit vec- 
tor in the momentum space. The crystallographic c-axis 
in the left half-space is selected parallel to the partition 
between the superconductors (along the z-axis in Fig.l). 
To illustrate the results obtained by computing the for- 
mula H2L)f) . we plot the current-phase diagrams for two 
different geometries. These geometries correspond to the 
different orientations of the crystals in the right and left 
sides of the interface (Fig.l): 

(i) The basal a6-plane in the right side has been rotated 
around the c-axis by a; 6i||c2. 

(ii) The c-axis in the right side has been rotated around 
the fr-axis by a; t>i||b 2 . 

Further calculations require a certain model of the gap 
vector (order parameter) d. 

IV. ANALYSIS OF NUMERICAL RESULTS 

In the present paper, the nonunitary /—wave gap vec- 
tor in the B— phase (low temperature T and low field H) 
of superconductivity in UPtj, compound has been consid- 
ered. This nonunitary bipolar state which explains the 
weak spin-orbit coupling in U Pt$ is |5j : 

d(T, v F ) = A (T)fc z (x (k 2 x - k 2 y ) + y2ik x k y ). (21) 
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FIG. 5: Current tangential to the interface versus the phase 
difference <f> for the junction between nonunitary bipolar 
/—wave superconducting bulks, T/T c = 0.15, geometry (ii) 
and the different misorientations (^component) . 



The coordinate axes x, y,z are selected along the 
crystallographic axes a, b,c in the left side of Fig^ 
The function Ao =Ao (T) describes the dependence of 
the gap vector on the temperature T (our numerical 
calculations are done at the low value of temperature 
T/T c = 0.1). Using this model of the order parameter 
(|21|l and solution to the Eilcnberger equations 12011 . 
we have calculated the current density at the interface 
numerically. These numerical results are listed below: 

1) The nonunitary property of Green's matrix diagonal 
term consists of two parts. The explicit part which 
is in the B mathematical expression in Ea. (|20|l and 
the implicit part in the f2i j2 and di j2 terms. These 
O12 and di :2 terms are different from their unitary 
counterparts. In the mathematical expression for Jli j2 
the nonunitary mathematical terms Ai j2 are presented. 
The explicit part will be present only in the presence of 
misorientation between gap vectors, 

B = id\ x d 2 • (Ai+A 2 )(?7fl 2 - e)(r}Q.\ + e), but the 
implicit part will be present always. So, in the absence 
of misorientation (di||d 2 ), although the implicit part of 
nonunitary exists the explicit part is absent. This means 
that in the absence of misorientation, current-phase 
diagrams for planar unitary and nonunitary bipolar 
systems are the same but the maximum values are 
slightly different. 

2) A component of current parallel to the interface j z for 
geometry (i) is zero similar to the unitary case [8( while 
the other parallel component j x has a finite value (see 
Fig^J. This latter case is a difference between unitary 
and nonunitary cases. Because in the junction between 
unitary /—wave superconducting bulks all parallel 
components of the current (j x and j z ) for geometry (i) 
are absent 0. 
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FIG. 6: Current tangential to the interface versus the phase 
difference (f> for the junction between nonunitary bipolar 
/—wave superconducting bulks, T/T c — 0.15, geometry (ii) 
and the different misorientations (z— component). 

3) In Figs l2l3l the Josephson current j y is plotted for 
a certain nonunitary model of /—wave and different 
geometries. Figs l2l3l are plotted for the geometry (i) and 
geometry (ii) respectively. They are completely unusual 
and totally different from their unitary counterparts 
which have been obtained in 

4) In Fig[21 for geometry (i), it is observed that by 
increasing the misorientation, some small oscillations 
appear in the current-phase diagrams as a result of the 
non-unitary property of the order parameter, . Also, the 
Josphson current at the zero external phase difference 
4> = is not zero but has a finite value. The Josephson 
current will be zero at the some finite values of the phase 
difference. 

5) In Fig|3 for geometry (ii), it is observed that by 
increasing the misorientation the new zeros in current- 
phase diagrams appear and the maximum value of the 



current will be changed non-monotonically. In spite of 
the FijJ21for geometry (i), the Josephson currents at the 
phase differences <j) — 0, (f> — 7r, and cf> — 2w are zero 
exactly. 

6) The current-phase diagram for geometry (i) and 
x— component (Fig0J) is totally unusual. By increasing 
the misorientation, the maximum value of the current 
increases. The components of current parallel to the 
interface for geometry (ii) are plotted in Fig[S] and 
Fij^JJ All the terms at the phase differences = 0, 
4> = 7T, and 4> = 2ir are zero. The maximum value of 
the current-pase diagrams is not a monotonic function 
of the misorientation. 

V. CONCLUSIONS 

Thus, we have theoretically studied the supercur- 
rents in the ballistic Josephson junction in the model of 
an ideal transparent interface between two misoriented 
UPt^ crystals with nonunitary bipolar /—wave super- 
conducting bulks which are subject to a phase difference 
</>. Our analysis has shown that misorientation between 
the gap vectors creates a current parallel to the inter- 
face and different misorientations between gap vectors 
influence the spontaneous parallel and normal Joseph- 
son currents. These have been shown for the currents in 
the point contact between two bulks of unitary axial and 
planar /-wave superconductor in Q separately. Also, it 
is shown that the misorientation of the superconductors 
leads to a spontaneous phase difference that corresponds 
to the zero Josephson current and to the minimum of 
the weak link energy in the presence of the finite sponta- 
neous current. This phase difference depends on the mis- 
orientation angle. The tangential spontaneous current is 
not generally equal to zero in the absence of the Joseph- 
son current. The difference between unitary planar and 
nonunitary bipolar states can be used to distinguish be- 
tween them. This experiment can be used to test the 
pairing symmetry and recognize the different phases of 
UPt 3 . 
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